Solving bifurcation diagrams using fixed points as control parameters 



David Andrieux 

We propose to determine the bifurcation diagrams of fixed points using their coordinates as control 
parameters. This method can lead to exact solutions to otherwise intractable bifurcation problems. 



EXECUTIVE SUMMARY 



We show that bifurcations diagrams can be obtained by using the dynamical variables as control parameters. Instead 
of looking at the evolution of a fixed point as a function of the system parameters, we deduce the system parameters 
as a function of the fixed point coordinates. 

This approach has several advantages. In many cases analytical solutions can be obtained as it bypasses the need to 
solve a nonlinear system. Furthermore, by construction, it provides a global analysis. The whole bifurcation diagram 
is obtained at one go. 

The key insight behind this result is that the system parameters are usually coupled linearly to the dynamics. 
Therefore, given a fixed point, the bifurcation problem is reduced to solving a set of linear equations for the system 
parameters. The whole bifurcation diagram is then constructed by varying the fixed point coordinates. 

We introduce the technique through examples. A similar reasoning can then be applied to other dynamical systems 
of interest. 

I. EXAMPLE 1. FITZHUGH-NAGUMO MODEL 

Our first example is the FitzHugh-Nagumo model 

— = v — v —w + 1 (la) 
t— — = v — aw . (lb) 

The FitzHugh-Nagumo describes a prototype of an excitable system (e.g., a neuron). 

We use this simple example to illustrate how to simplify the bifurcation analysis by avoiding solving a cubic 
nonlinearity for each values of the parameters / and a. 

The key observation is that / and a are coupled linearly to the dynamical variables v and w. Therefore, a fixed 
point (v, w) will impose a set of linear constraints for the parameters (J, a). Solving this linear system leads to 

I(v,w) = w — v + v 3 (2a) 
a(v,w) — v/w. (2b) 

These equations give the combination (I, a) for which (v,w) is a fixed point of ([I]). We can now obtain the whole 
bifurcation diagram analytically by varying the fixed point (v, w) and calculating its stability. Notably, this is achieved 
by solving a linear system rather a cubic nonlinearity as in the standard approach. 

II. EXAMPLE 2. BISTABILITY IN POPULATION NETWORKS 

We now consider a system of two variables x\ t2 G [0, 1] interacting through a positive, increasing function F: 

ax i x i 

— - = -+ 1 (l-x 1 )F{W 11 x 1 + W 12 x 2 ) (3a) 

at t 

^ = -—+j(l-x 2 )F{W 21 x 1 + W 22 x 2 ), (3b) 
at t 

where the coupling coefficients satisfy 

Wn = W 22 and W 12 = W 21 . (4) 
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Such dynamical systems appear, for example, in mean-field descriptions of neural networks (see, e.g., RefP). 

This system admits at least one symmetric fixed point (xo,Xo). Depending on the coupling strengths, the network 
can also admit asymmetric states. However, none of the fixed points (symmetric or asymmetric) can be obtained 
analytically. 

We use our method to study the bistability properties of the network (1 3 h in terms of the coefficients Wm. The 



!ior< 



interaction between the dynamical variables x% and the Ws is linear. Therefore we can determine the Ws compatible 
with a fixed point (xi,x 2 ). Noting that F is invertible since it is increasing, we have that 



xi x 2 \ ( Wii 
x 2 x\)\ W 12 



K 2 



where 



The solution of this system reads 



Ki = F~ 



7t(1 - xCj 



Wn(xi,x 2 ) = W 22 (xi,x 2 ) 
W 12 (xi,x 2 ) = W 2 i(xi,x 2 ) 



XyKy - X 2 K 2 

(xi + x 2 )(xi - x 2 ) 

X\K 2 - x 2 K\ 
(xi +x 2 )(xi - x 2 ) 



(5a) 
(5b) 



when x\ ^ x 2 . 

Equations ([5| express the coupling parameters in terms of the fixed point coordinates (xi,x 2 ). We can now obtain 
the whole bistability region analytically by varying these coordinates. This bistability region is depicted in Fig[T] 
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FIG. 1: Bistability region of the population network ([3]). The colored region depicts the couplings for which symmetric 
and asymmetric stable fixed points coexist. The color code indicates the value maxjii , X2) of the asymmetric state. The 
function F as well as the values of 7 and r used in this figure can be found in RefE The coupling coefficients are determined 
by Eqs. j5j. 



Note that, in this example, the standard problem of finding the fixed points as a function of the coupling is not 
analytically tractable. Nonetheless, our approach allows us to obtain the complex bifurcation diagram seen in Fig[l] 
exactly. Furthermore, it gives the global solution, i.e. the whole diagram, at one go. In particular, there is no need 
to follow different fixed points or branches to reconstruct the diagram. 
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